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Abstract 

We theoretically investigate the use of Bessel light beams generated using 
axicons for creating optical dipole traps for cold atoms and atomic waveg- 
uiding. Zeroth-order Bessel beams can be used to produce highly elongated 
dipole traps allowing for the study of one-dimensional trapped gases and re- 
alization of a Tonks gas of impentrable bosons. First-order Bessel beams are 
shown to be able to produce tight confined atomic waveguides over centimeter 
distances. 
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I. INTRODUCTION 



The experimental realization of Bose-Einstein condensation (BEC) in dilute alkali vapors 
0II has generated substantial theoretical and experimental research activity ||. BEC 
provides an important quantum system where one can observe various phenomena such as 
superfluidity ||], vortices |5]||, spin domains @ or soliton-like behaviour |^-[T0|. Central to 
future advances with these quantum degenerate systems is their manipulation by external 
potentials. The use of light fields offers an exciting avenue in this respect as optical potentials 
can be state independent, offer an array of spatial forms and can be rapidly switched. They 
offer very good prospects for novel dipole traps and matter-wave guides. 

Previous work has been reported for optical dipole traps based on standard Gaussian 
light beams JTT]-|T^|. Free-space propagating light beams such as Laguerre-Gaussian light 
beams |14[ and Bessel light beams |H| are excellent candidates for advanced all-optical 
manipulation of quantum gases. Laguerre-Gaussian light beams offer a dark hollow central 
region for guiding and focusing atoms f[6|,|PT|| . Their annular form can also be used to 



realise toroidal optical dipole traps for BEC [fLql , enabling studies of persistent currents. 
Bessel light beams are solutions of the scalar Helmholtz equation that are propagation 
invariant. This immunity to diffraction coupled with the small size of their central region 
means they offer unique characteristics in the optical domain. In this work we discuss the use 
of Bessel light beams for trapping and guiding cold atoms at or close to quantum degeneracy. 
Specifically we study the Bessel light beam optical dipole trap with the aim of generating one- 
dimensional quantum gases. In particular we look at the possible realisation of a Tonks gas of 
impenetrable bosons which exhibit fermionic like excitations jn|. Furthermore, we analyse 
in detail the waveguiding of a matter-wave beam along a higher-order Bessel light beam. 
This method could provide an important route to realising all-optical atom interferometers 
and further the loading of magnetic waveguides. 



II. BESSEL LIGHT BEAMS 



In this Section we review the basics of the generation of Bessel light beams. Ideal 
forms of such beams are impossible to realise as they would have infinite extent and carry 
infinite power 
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However, finite approximations to these beams can be realised that 
propagate over extended distances in a diffraction free manner. Holographic methods offer 
efficient generation of Bessel beams |20|j2l[1 . Another very efficient method to generate an 
approximation to a Bessel beam is by use of a conically shaped element termed an axicon and 
it is this form of Bessel beam generation we concentrate upon in this work. This generation 



method is discussed in detail elsewhere p2"||2"3|j but we review it here for completeness and 
clarity in notation. Our main interest is using zeroth-order Bessel beams for optical dipole 
traps and first-order Bessel beams for atomic waveguiding, and we shall concentrate on these 
cases. 

Bessel beams are solutions of the free-space wave equation which propagate with un- 
changing beam profile along the propagation axis which we take as z in cylindrical coordi- 
nates (r, 8, z). The electric field of a monochromatic, linearly polarized ideal Bessel beam of 
order t and frequency u>l is [ pT5| , pOH 22 
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E(r, 9, z,t) = ~ (E J l (k r r)e« k ** +iB - u rt + c.c.) , (1) 

where Eq is a scale electric field value, Jt is the £th order Bessel function, i > is the 
azimuthal mode number which we take as positive for simplicity in notation, and k r and k z 
are the radial and longitudinal wavevectors such that k 2 = k 2 + k 2 z with k = ujl/c = 2tt/\l. 
The zeroth-order solution Jo has a central maximum surrounded by concentric rings of 
roughly equal power while the higher-order solutions Ji have zero on-axis intensity also with 
concentric rings. 

The method of generating a Bessel beam by use of an axicon is based on the observation 
that the Fourier transform of the Bessel beam solution ([!]) over the transverse plane (r, 9) 
is an infinitely high ring in the spatial frequency domain peaked at K = k r , where the 
phase varies from zero to 2tc£ around the peak: the Bessel beams may therefore be viewed 
as a superposition of plane-waves with transverse wavevectors lying on a ring of magnitude 
k r . A finite realization of such a beam may be produced by passing a Laguerre-Gaussian 
(LG) beam of order i through an axicon as illustrated in Fig. [I]. The axicon is a conically 
shaped optical element which imparts a phase shift ax (r, 9) = k ax r to an incident field where 
k ax = k(n— 1)7, n being the refractive index of the axicon material, and 7 the internal angle 
of the element. By choosing k ax = k r for a specific Bessel beam the axicon imposes the 
ring of transverse wavevectors characteristic of the Bessel beams on the incident beam, and 
these plane-waves come together past the axicon to produce a Bessel beam. Single-ringed 
LG beams (i.e. with radial mode index p = 0), which may be produced using holographic 
elements, have electric field envelopes (plane-wave factor exp(i(kz —u)z,t)) removed) at focus 



of the form 14 



2P /2r^ 



= °> = V;^U eX P^>o)e^, (2) 



where Po is the power at z = incident on the axicon, wq is the Gaussian spot size, and i is 
the azimuthal mode number. Then a LG mode of order i incident on a axicon with k ax = k r 
has the appropriate azimuthal phase variation to produce a £th order Bessel beam past the 
axicon. In particular, from geometrical considerations and for an LG mode of size Wq, we 
expect the ring of plane-waves imposed by the axicon to overlap spatially over a longitudinal 
range z max sin(i9) w w , with sin(i9) = k r /k. This gives the estimate [|23| 



kw 



r 



for the range over which the plane-waves overlap and produce a finite realization of the 
Bessel beam profile Ji(k r r) past the axicon: The larger the input spot size w the larger 
the center portion of the transverse plane over which the actual field approaches the ideal 
Bessel beam. 

Numerical simulations and accompanying experiments have verified this physical picture 
for Bessel beam generation using an axicon p2| , |23|| . The numerical simulations used the 



Fresnel diffraction integral to propagate the input LG field (eq. @) times the phase aber- 
ration exp(icj) ax (r,6)) due to the axicon, to distances beyond the axicon. Figure |2| shows a 
gray-scale plot of the field intensity I(x, y — 0, z) as an example of this numerical propaga- 
tion for I = 1,Wq — 0.29 mm, = 780 nm, n = 1.5, and 7 = 1°, giving k r = 6.8 x 10 4 
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m _1 . For this example z max = 3.4 cm and the simulation shows that this indeed estimates 
the scale of the Bessel beam propagation. As expected for a J\ beam the intensity profile 
has a dark fringe at the center of width of about wb = l/k r = 14.7 /iin. 

The evolution of the field past the axicon may also be approximated using the method 
of stationary phase applied to the Fresnel integral |23| , which yields the following expression 
for the field intensity h(r, z) for an £th order LG input mode 

h(r,z) « — M U- ) exp(-2z 2 /^axU 2 (M, (4) 

with 2 max given by Eq. (|3|). Figure |3| shows the on-axis intensity 1(0, z) versus propagation 
distance z past the axicon obtained using both the Fresnel integral approach (solid line) and 
the approximate solution (dashed line) for the same parameters as Fig. ^] except t = 0, 
with excellent agreement (for 1 = the intensity has a peak on-axis). In particular, by 
differentiating Eq. @ we find that the peak intensity evaluated over the whole transverse 
plane occurs at 



_ V21+T 

-^peak 2 ^max) \p) 

which gives z pea k = 1.7 cm, in excellent agreement with the results in Fig. |3|. The approxi- 
mate expression (H) is applicable over the center region of the Bessel beam if z > kjk 2 r p3 |. 
Using Eq. (|3|) this condition can be written as z > 2 max / (kw^): Recognizing that kw\ = 2zr 
is twice the Rayleigh range zr of the input LG beam |l4|, and = l/^ r is a measure of 
the width of the central lobe of the Jo Bessel beam or the central dark fringe for a J\ beam, 
we see that if 

Zma * 2WB « 1, (6) 



Zr w 

then Eq. (|j) should be valid over most of the propagation range of the Bessel beam around 
the position of the peak z pea ^. That is, under conditions where the central spot or dark fringe 
of the Bessel beam is narrow compared to the input LG spot size, which is what we want, the 
stationary phase approximation should be valid. For the above examples wb/wq = 1/19.7, 
and this is typical of parameters we consider. We shall assume this condition is satisfied 
hereafter and use the approximate expression (|J) in the remainder of this paper. 

Finally, we point out one further feature of Bessel beams that highlights their utility, 
namely, that their radial width wb = l/k r = l/(k(n — 1)7) is determined solely by the laser 
wavelength and the axicon parameters, whereas their longitudinal extent z max = kwo/k r = 
kwoWs is also dependent on the incident LG spot size Wq. This means that the longitudinal 
extent and radial confinement can be varied independently. This is in contrast to a Gaussian 
beam of equal spot Wo for which the longitudinal extent of the focus is the Rayleigh range 
zr = kwl/2. 



III. BESSEL OPTICAL DIPOLE TRAPS 

In this Section we examine the use of zeroth-order or J Bessel beams for creating elon- 
gated optical dipole traps for BECs, and assess their utility for realizing one-dimensional 
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trapped gases and a Tonks gas [19[] of impenetrable bosons f24] , |25 |. Here we concentrate on 
the case of bosonic atoms but similar considerations apply to degenerate fermionic atoms 
| 26| in Bessel beam traps. 



A. Gross-Pitaevskii equation 

The Gross-Pitaevskii equation (GPE) describing the macroscopic wave function ip(r,t) 



of an iV-atom BEC in an optical dipole trap can be written as [27 



dib h 2 

ih ^t = ~2M V ^ + y(r ' + U ° N ^> ( ? ) 

where V 2 = V 2 + d 2 /dz 2 is Laplacian which is the sum of the radial and longitudinal 
Laplacians, M is the atomic mass, U$ = AnTi a/M is the effective three-dimensional short- 
range interaction strength, with a being the s-wave scattering length. The potential term 
on the right-hand-side 

= (8) 



8A 



'Sat 



describes the optical dipole potential with A = ujl — uja the laser detuning from the optical 
transition frequency u>a, T the natural linewidth of the optical transition, 7g at is the resonant 
saturation intensity, and I(r,t) = |e c|-E(r, z)\ 2 . For a red-detuned laser the potential is 
negative and the atoms are attracted to the regions of high intensity, whereas for a blue- 
detuned laser the atoms are repelled into the low field regions. 



B. Optical dipole potential 

Here we investigate the properties of an optical dipole trap formed using a red-detuned 
Jo beam, so that the atoms are attracted to the intense central maximum of the intensity 
distribution given by Eq. (|J) with t = 

J (r, z) « 2**^0 f-TAjl (— ) exp(-2* 2 /4J J 2 (k r r). (9) 

We are primarily interested in tight bound optical dipole traps so we shall approximate the 
optical dipole potential using the parabolic approximation to the full intensity Io(r,z) 
around the intensity maximum at (r = 0, z = ^ pea k): 

V(r, z) - V(0, z peak ) w l -Mn 2 r0 (r 2 + \ 2 {z - ^ pcak ) 2 ) , (10) 

with 

, , s hT 2 P k , „ , 2y/2 w B , \ 

n 2 r0 = exp -1/2 — y— y k 2 , A = — ^— = 2.83—^. 11 

r0 PV 1 4 A M/sat 2 max r ' kw z max 
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A red-detuned (A < 0) Jo optical dipole potential therefore provides confinement in both 
the radial and longitudinal directions. Here Q r o is the radial oscillation frequency with cor- 



responding ground state oscillator width w r o = yh/MQ r o, and A is the ratio between the 
longitudinal and radial trap frequencies Q z o/Q r o = A [^-|3(J, which also determines the 



aspect ratio between the radial and longitudinal ground state widths w t q/w z q = y/X (in the 
absence of many-body repulsion). We remark that even for tight transverse confinement an- 
harmonic corrections beyond the parabolic approximation for the longitudinal (z) variation 
of the optical dipole potential ([II]) are generally required. This is evident from Fig. |3] where 
the on-axis intensity variation for the Jo Bessel beam is not symmetric around the peak. For 
this work we shall restrict ourselves to the parabolic approximation. We remark, however, 
that anharmonic corrections to the parabolic approximation can have pronounced effects 
on the trap properties such as the condensate fraction and the frequencies of the collective 
excitations that can be supported by Bessel atomic waveguides, especially for shallow traps 

So far we have ignored the effects of gravity which we assume acts in the radial direction 
so the longitudinal axis of the trap is horizontal. For tight radial confinement gravity will 
serve simply to displace the origin of the radial motion, and this applies under conditions 
such that Mil 2 w 2 /2 ^> Mgw r0 with g = 9.81 ms" 1 . Using w r0 = J%/MQ r o, this lead to 



the condition on the confinement of the ground state, io r o ^ y h 2 /(2M 2 g). For example, 
for 85 Rb the radial confinement has to be less than 0.31 /mi, whereas for 23 Na a confinement 
w r0 -C 0.74 /im would be sufficient. We assume that this condition is satisfied and hereafter 
neglect the effects of gravity. 



C. One-dimensional trapped gases 

From Eq. (|TT| ) the parameter A is seen to determine the anisotropy of the Jo optical 
dipole trap. It is well known that in the limit A 1 highly asymmetric cigar-shaped traps 



are formed which are elongated along the longitudinal direction pq-|30|. Bessel beams are 
exceptional in this regard as they can produce extreme asymmetries. To illustrate this 
recall that wo is the input LG spot size before the axicon, for example wq = 1 mm, so that 
A = 4.8 x 10~ 4 for a A^ = 1064 nm. We stress that one obtains such large asymmetries 
using Bessel beams without sacrificing the radial confinement, as the radial (wg = l/k r ) and 
longitudinal (z max ) extents of the Bessel beam are independently variable. In contrast, for a 
Gaussian optical dipole trap changing the focused spot size wq also changes the longitudinal 
extension of the dipole trap set by the Rayleigh range zr = kw\j1 |H 



The limit A <§; 1 corresponds to the regime of one-dimensional trapped gases pq , |29|j3C 
in which the radial variation of the macroscopic wave function is effectively frozen as the 
normalized ground state mode u g (r) of the radial optical dipole potential 

E g u a = -^V r Y, + ^MQ 2 r0 r\. (12) 

Writing the macroscopic wave function in a form reflecting the single-radial mode nature of 
the solution 

r(>(r,z,t)=u g (r)<f>(z,t)e- iE ° t '\ (13) 
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and combining Eqs. fl7|), fllCf) , (|12"D , and (p~3|) we obtain the one- dimensional Gross-Pitaevskii 
equation for the reduced system |2^]30| 



(14) 



where 



g = Uo J 2nrdr\u g (r) 
2hQ r Qa, 



2h 2 a 
Mw 2 



(15) 



rO 



is the effective one-dimensional short-range interaction strength p5| . In steady-state we set 
(p(z, t) = x( z ) exp(—ifit/h) in Eq. ([14]) with \i the chemical potential of the one-dimensional 
system. Then using the Thomas-Fermi approximation pS| for high density in which the 
kinetic energy term is neglected p8] |, Eq. ([14]) yields for the one-dimensional density 



p(z) = \ X (z)\ 2 = ^ 1 



( Z ^peak ) 



| Z ^peak | ^ %r, 



(16) 



with 



/3#iV\ 2/3 fMtt 2 r0 \ 



2\ V3 



V 4 



') 



3Nw 2 a 
A 2 



1/3 



(17) 



This density solution has a peak one-dimensional density p pea k = 1.5(N/L Z ) oc N 2 ^ 3 with a 
longitudinal length L z = 2z m . 

The Gross-Pitaevskii equation ([14]) has previously been investigated as a model for a one- 
dimensional Bose-Einstein condensate in a number of situations, including the ground state 
32[ and dynamics of cigar shaped traps [|34] -|36]j, dark solitons [p0] ,p|, p7l , |38| , bright solitons 
for negative scattering lengths |]29| , |30| , |39[ | , gap solitons in optical lattices |§, atom waveguides 
and as a model Luttinger liquid [13]. Here our goal is to highlight the utility of J 



optical dipole traps for experimental studies of one-dimensional trapped gases. To illustrate 
the basic scales involved in these elongated traps we consider the case of the = 780 nm 
transition of Rb with /sat — 16W/m 2 and r = 2tt x 6.1 MHz. A far red-detuned laser at a 
wavelength of Al = 1064 nm is used to generate a Bessel beam with a longitudinal extent 
2 max = 10 cm and a central spot of full- width radius 3 /im (corresponding to wb = 1-25 /zm 
and consequently u7o = 13.6 mm). If a laser power of Pq = 5W is used this results in a 
trap potential of about — 49/iK depth with radial trap frequency Q r o = 2tt x 8.8 kHz, radial 
confinement of w r o = 82 nm, and an aspect ratio of A = 3.5 x 10~ 4 . Then for iV = 10 4 
atoms and a scattering length a = 5 nm, for example, we find L z = 2z m = 1.9 mm, and a 
peak density p pca k = 6.4 x 10 4 cm -1 . For comparison, a focused Gaussian with the same 
radial confinement (and trap depth) gives a trap with an aspect ratio of only 8 x 10 -2 . We 
remark, however, that the power required to obtain the same radial confinement for the 
Gaussian beam is much less, typically on the order of milliwatts |TTJ. This arises because 



for the Bessel beam the trapped atoms only experience the intensity in the central peak of 
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the beam, and the concentric rings surrounding the peak, though key in realizing the long 
extent of the Bessel optical dipole trap, do not directly affect the atomic confinement. 

Bessel beam traps will suffer from the same loss mechanisms as other shallow far-off 
resonance dipole traps. Apart from heating due to diffractive collisions with background 



gas the trap lifetime will be limited by laser-noise-induced heating |4l! |. The expected 
loss rates due to intensity fluctuations and pointing instability of the laser beam for a Bessel 
beam trap should be comparable to those in Gaussian beam traps with identical radial 
confinement and trap depth. 

Finally, we note that the Jo optical dipole trap provides a means to accelerate an atomic 
trapped one-dimensional gas. Consider a one-dimensional gas in the ground state of the Jq 
trap peaked around z = z pca k = kwo/2k r . If we now start to slowly increase the LG spot 
size wo(t), then the longitudinal position of the peak density of the trapped gas should vary 
in time according to z vea _^{t) = kw Q (t)/2k r . Such a scheme provides a means to impart a 
longitudinal velocity to an initially stationary trapped gas. 



D. Realization of a Tonks gas 

In the above discussion we tacitly assumed that the system of cold atoms formed a 
BEC. However, Petrov et al. have theoretically studied the diagram of state for a one- 
dimensional gas of trapped bosons, assuming A < 1, and found that a true BEC, or at 
least a quasi-condensate, with concomitant macroscopic occupation of a single state, is only 
attained for high enough particle numbers N > with 



For N < and temperatures T < NhQ z0 , one obtains a Tonks gas of impenetrable bosons 
for which hard core repulsion between the bosonic atoms prevents them from penetrating 
through each other in the one- dimensional system, and the system acquires properties rem- 
iniscent of a one-dimensional system of fermionic atoms. This remarkable property of the 
Tonks gas is related to the breakdown of the spin-statistics theorem in one-dimension, and is 



reflected in the Fermi-Bose mapping for this system first elucidated by Girardeau ||46fl , and 
applied to atomic waveguides by Olshanii P3. Recent theoretical investigations of Tonks 
gases have shown that they can support dark soliton structures JlTlKBL and also that their 



coherence properties are significantly different from the corresponding BEC flSJSTl. Fur- 



thermore, recent experimental developments suggest that Tonks gases should be realizable 
in magnetic atom waveguides |40|-|42|, and Bongs et al. |36|] have proposed a hybrid trap 



composed of optical dipole trap formed with a first-order LG beam combined with magnetic 
longitudinal trapping. Here we examine the utilty of Jo optical dipole traps for realizing a 
Tonks gas. 

The highly elongated Bessel beam discussed in the previous section (A^ = 1064 nm, 
Po = 5W, z max = 10 cm, wb = 1.25 /mi) would be an ideal candidate for the experimental 
realization of a Tonks gas. The low aspect ratio A = 3.5 x 10~ 4 and tight radial confinement 
w r o = 82 nm result in a high upper boundary iV* for the particle number of the Tonks gas. 
For the commonly used 87 Rb isotope with a scattering length a = 5 nm one finds = 420. 
Although this is still a fairly low value it should be possible to experimentally realize a small 
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87 Rb Tonks gas. However, more promising would be the use of the 85 Rb isotope, where a 
Feshbach resonance can be used to tune the normally negative scattering length to positive 



values of several hundred nanometers magnitude ||5l| . As iV* is proportional to the square 
of the scattering length even a moderate increase to a = 50 nm would make it possible to 
create a larger Tonks gas, with say N = 2000 atoms, which should be easily detectable. 

The Bessel beam trap offers some advantages compared to alternative suggested ap- 
proaches using magnetic waveguides |]4"0"|-fl2]] and a hybrid magnetic-optical trap Firstly, 
it involves only a very simple all-optical system for which the aspect ratio of the trap may 
be controlled simply by varying the Gaussian spot size incident on the axicon. More specif- 
ically, being an all-optical, it does not involve material surfaces as in magnetic waveguides, 
which can cause matter- wave decoherence p2| , |53| . Furthermore, it allows for the possibility 
of trapping multiple magnetic sublevels and the investigation of multi-component Tonks 
gases, which would not be possible in the hyrbrid magnetic-optical trap of Bongs et al. p6 



IV. BESSEL BEAM ATOMIC WAVEGUIDES 

The higher-order Bessel beams Je with £ > have zero on-axis intensity surrounded by 
intense concentric rings, and a blue-detuned beam can trap therefore atoms radially in the 
dark hollow core of the beam. In addition, since the intensity vanishes on-axis the higher- 
order beams produce negligible longitudinal confinement (see Fig. ^| for the J\ Bessel beam) 
in comparison to the Jo Bessel beam. The higher-order Bessel beams therefore present 
the opportunity to realize atomic waveguides with tight radial confinement over distances 
comparable with the propagation range z max of the Bessel beams. Furthermore, the J\ beam 
is special in that it provides a parabolic radial trapping potential to lowest order (see below), 
whereas the higher-order trapping potential vary as r 2e near the axis. We shall therefore 
concentrate on the Ji atomic waveguide. 



A. Atomic waveguide potential 

Here we investigate the properties of an atomic waveguide formed using a blue-detuned 
(A > 0) Ji beam, so that the atoms are repelled into the zero intensity central minimum of 
the intensity distribution given by Eq. (|J) with I = 1 

h(r, z) » AtxKw ( —Y e-^l*™ J*(k r r). (19) 

Substituting this expression in the optical dipole potential @ we write the atomic waveg- 
uiding potential in the form 

V(r, z) = \Mtf rl {z) (4^|^) - \Mtf rl {z)r\ (20) 

where in the last form of the potential we used the approximation J\{x) ~ x/2 which is 
applicable for tight radial confinement. Here the ^-dependent radial oscillation frequency 
Q r i(z) is given by 
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/ \ 3/2 

nn(z) = n r i{z pcak ) x — exp(-0.75(z 2 /^ cak - 1)), (21) 

\ %>eak / 

with Zpeak — V3z mSiX /2 the position of the peak intensity for the J\ Bessel beam, and 
^ri(-2 P oak) is the peak radial oscillation frequency given by 

tft(w)=(^-J ^(-3/2)^^— e (22, 
In comparison to the radial oscillation frequency of the Jo Bessel beam trap we find 

^rl(^pcak) ~ 0.5Q r o. 

The Ji Bessel beam therefore defines an atomic waveguide whose radial confinement 
peaks at z — z pca k and varies with the longitudinal coordinate z. If we take the effective 
length L of the waveguide to be the full width at half maximum of O r i(z) versus z then by 
inspection of Eq. (^T|) we find L m z max , as may have been anticipated physically based on 
the properties of the Bessel beams. 



B. Gross-Pitaevskii equation 

To examine the properties of a BEC propagating in a J\ atom waveguide we shall ap- 
proximate the waveguide as invariant along the z-axis over the length L = 2 max and calculate 
the ground state radial mode of this system. In general we should solve the propagation 
problem of the matter wave field through the varying atomic waveguide potential, and how 
the atoms are funneled into the waveguide, but that shall be the subject of future a paper. 

To proceed we consider the GPE ([/]) for a BEC of momentum p z per atom moving along 
the z-axis in a cylindically symmetric J\ Bessel atom waveguide which is invariant along the 
z axis. Then writing the macroscopic wave function as 



N^(r, z, t) = y?(r)e i{p ^- M)/? \ (23) 




/* - ^7 I V = ~Wn[lhi + ZlL)<P + VoJf(k r r) V + f/oM V (24) 



the GPE becomes 

2M 

where \i is the chemical potential, and Vq = 2M^ 1 (z pea j ! )/^ characterizes the strength of 
the Ji waveguide. The wave function over the atomic waveguide effective length L = z mSuX 
is normalized to the number of particles 

z max / 2Tirdr\ V {r)\ 2 = N. (25) 
Jo 

To facilitate numerical solution of the GPE we introduce scaled variables as follows |)4 

C = Kr, <p(r) = J^MO, (26) 
in terms of which the GPE for <p(C) becomes 
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a( P = - 375 + 7T77 ^ + (OV + c^l V (27) 



Here £ r = h 2 k 2 /2M sets the energy scale, a = (/i — p 2 /2M)/E r is the scaled energy 
eigenvalue, (3 = V /E r , and c = ±1 sets the sign of the many-body interactions. The wave 
function <p(C) is now normalized as 



o 



CO 

(d(\ V (()\ 2 = V N = n, (28) 



where the dimensionless variable rj = (k 2 \U \/ L)/(2nE r ) = A\a\/L. 

Strictly speaking the J\ atomic waveguide potential in the GPE fl2"T| ) does not have bound 
state solutions, since any initial wave function localized in the central minimum of the J\ 
potential will ultimately tunnel out over the finite potential barrier due to the first peak 
of the Bessel beam. However, for tight confined atoms the tunneling time can be made 
arbitrarily long, and here we ignore tunneling to lowest order. Figures ^ and show the 
results of the numerical solution of the GPE (^) for (3 = 10 3 and c = 1, that is repulsive 
many-body interactions. These numbers would for example describe a guide for rubidium 
with a radius of 5 fim (corresponding to wb = 2.78 /im) and length L = 5 cm using 40 mW 
of light at Xl = 1064 nm. The numerical method used to solve eq. (^7]) was the same as in 
Ref. [53|. Figure ^ shows the variation of the scaled energy eigenvalue a versus the scaled 



number of particles n = r/N, and we see that as n increases so does the energy due to the 
repulsive many-body interactions. For n — > the scaled energy reduces to that of the atomic 
waveguide potential (see below). Figure § shows the scaled wave function ip(() versus scaled 
radial coordinate ( for n = 0.88 (a = 50) (solid line), n = 13.4 (a = 150) (dotted line), and 
n = 43 (a = 250) (dashed line). Here we see that as the scaled number of atoms increases so 
does the width of the wave function, as expected physically. Clearly, there is a limit to the 
allowable width of the wave function, and hence the number of atoms, as the tunneling out 
of the atomic waveguide alluded to above will become more relevant as the wave function 
width approaches that of the central minimum of the J\ Bessel beam. 



C. Thomas-Fermi approximation 

In this section we discuss approximate solutions to the GPE (|2~7D to provide a framework 
for the numerical solutions. First for tight confinement, so that tunneling out of the trap is 
negligible, we approximate Ji(C) ~ C/2, giving 

( d 2 1 d \ (3( 2 , l2 

For this to be valid the macroscopic wavefunction should not extend beyond the first peak 
of the J\ Bessel function at ( c = 1.8. Then in the limit of a small number of atoms r\ N — > 
the GPE has the Gaussian ground state solution 



^(C)ocexp(-C7Co), a = V& ( 30 ) 
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with Co = 2/3 x / 4 : this is only valid if Co < Co or /3 3> 1. For (3 = 10 3 this gives a = 31.6 in 
agreement with Fig. |] as rjN — > 0. In the limit of large 7]N we can use the Thomas- Fermi 



solution [£8] in which the kinetic energy term is neglected. This yields the solution 



\<p(0\ = vW i-CVC , C<U = 2 V > ( 31 ) 




and «tf = y/rjWp. This solution is only valid if ( m < ( c = 1.8, which for a given value of j3 
places an upper bound on the scaled number of atoms n 

n = r]N < (3. (32) 

In Figs. H and [| we have restricted n < 50, in accordance with the above upper bound, 
and the wave functions in Fig. [5] are all vanishingly small in the region ( > ( c = 1.8. As 
an example, for n = 43, /? = 10 3 the Thomas-Fermi solution predicts ( m = 0.9, which is 
smaller than the spatial extent of the wave function (dashed line) in Fig. [|, and q?tf — 207 
in comparison to a = 250 from the exact solution in Fig. ^. This discrepancy between the 
exact and approximate solutions is not surprising as we used the parabolic approximation 
to the optical potential in the Thomas-Fermi solution. However, the Thomas-Fermi theory 
captures the trends of the solution. 



V. SUMMARY AND CONCLUSIONS 

Bessel light beams have unusual properties in the optical domain. They have an immunity 
to diffraction over extended distances and offer an elongated and narrow central region. 
These features give them significant advantages over standard Gaussian light beams for 
the studies presented in this work. We have shown that zeroth-order Bessel light beams 
generated by use of an axicon offer an excellent method by which to generate optical dipole 
traps for one dimensional quantum gases. Typically, the ratio A between longitudinal and 
transversal trap frequencies for such a trap realised with a Bessel beam is nearly two orders of 
magnitude smaller than that that can be achieved with a Gaussian beam. Furthermore, we 
have shown that the Bessel beam offers a potential route for an experimental realisation of a 
Tonks gas of inpenetrable bosons. Additionally we have studied the waveguiding properties 
of a quantum gas along a J\ (hollow) Bessel light beam. The non- diffracting nature and small 
central minimum size here make it an excellent all-optical waveguide. Such waveguiding 
could be used to realise velocity filtering of cold atoms, all-optical atom interferometers and 
also offer a route to load magnetic waveguides with quantum degenerate samples. 
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1-0806, and the Department of Army Grant No. DAAD 19-00-1-0169. JS acknowledges 
support from an NSF graduate student grant. KD acknowledges the support of the UK 
Engineering and Physical Sciences Research Council and the Leverhulme Trust. 
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FIGURES 



FIG. 1. Illuminating an axicon with a LG mode of order t produces a Bessel beam of the same 
order within the shaded region. 

FIG. 2. Gray scale plot of the intensity I(x, y = 0, z) for an input LG beam with £ = 1 showing 
Bessel beam formation over an elongated distance. Other parameters are u>o = 290 ^m, Xl = 780 
nm, n = 1.5, and 7 = 1°, giving z max = 3.4 cm, and z pea k = 2.95 cm. 

FIG. 3. On-axis intensity variation for an input I = LG beam using both the exact Fresnel 
integral approach (solid line) and the expression (|j) based on the stationary phase approximation 
(dashed line). All other parameters are the same as Fig. 

FIG. 4. Scaled energy eigenvalue a versus scaled particle number n = rjN for (3 = 10 3 . Phys- 
ically, a is the chemical potential minus the longitudinal kinetic energy per particle scaled to the 
energy E r , and (3 = Vq/E t measures the depth of the Bessel optical potential scaled to E r . 

FIG. 5. Scaled transverse eigenmode </?(C) versus scaled radial coordinate £ for [5 = 10 3 and 
n = 0.88 (a = 50) (solid line), n = 13.4 (a = 150) (dotted line), and n = 43 (a = 250) (dashed 
line). 
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